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^ I I Abstract: In this paper we present a method to design paraunitary polyphase matrices of critically sampled rational 

■ filter banks. The method is based on (P, Q) shift-invariant systems, and so any kind of rational splitting of the frequency 

tyj I spectrum can be achieved using this method. Ideal (P, Q) shift-invariant system with smallest P and Q that map of a 

O . band of input spectrum to the output spectrum are obtained. A new set of filters is obtained that characterize a (P, Q) 

I shift-invariant system. Ideal frequency spectrum of these filters are obtained using ideal (P, Q) shift-invariant systems. 

T-H ■ Actual paraunitary polyphase matrices are then obtained by minimizing the stopband energies of these filters against 

^ ' the parameters of the paraunitary polyphase matrices. 

00 ; 
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' I. Introduction 

o : 

, A. Motivation 

• , Rational Filter Banks (RFB) are now well known tools for nonuniform subband decomposition of the signal. Given 

N rational numbers {^i • ' ' 7 ^"'^ }, where X^il^o^ ^ ^' '■^^ ^^^^ °f ™ ^ channel RFB is to split the frequency 
' spectrum [0, tt] of the signal into N nonuniform bands, where n-th band is given by E"=o f^"^' S"=!o^ f^^]- 

' An RFB uses a P/Q sample rate changer, shown in Fig.© in each channel. 

Such RFB are treated in detail in fTI and it is shown that all kind of rational band splitting can not be achieved 
by the RFB based on the above P/ Q sample rate changer, even if the ideal real coefficients filters are used. 

If P and Q are coprime, the structure shown in Fig.([T|) can be converted into a structure shown in Fig.©. In 
this figure Hp(z) is a P x Q polyphase matrix that consists of polyphase components of H{z). It is shown in |3] 
that this structure is more general than the previous structure, and if P and Q are allowed to have common factors, 
then any kind of rational splitting of the input frequency spectrum can be achieved. 

Such structures are called (P, Q) shift-invariant systems, since shifting the input by Q samples results in shifting 
the output by P samples. Many properties of (P, Q) shift-invariant systems are studied in H. It is shown that such 
structures can be represented by a P x Q LTI polyphase matrix Hp(z). In other words, let input x{z) and output 
y{z) be written in polyphase form as x{z) = X^iLo^ z^^Xi{z'^) and y{z) = X]i=o^ 2 *yi(^^) respectively, and let 



[xo{z),Xl{z), ■ ■ ■ ,XQ_l(2) 



Vpi^) = [ya{z),yi{z),- ■ ■ ,yp-i{z)]'^ 



(1) 



DRAFT 




Fig. 1. A typical brancli of a rational filter bank 




Fig. 2. Polyphase structure for RFB implementation 



then 

yp{z) = Hp(z)xp(z). (2) 

Let Sx{uJi,uj2) represent the frequency spectrum of the signal x{n) between frequencies uti < lu < uj2 ■ 
The {P, Q) shift-invariant systems described above can be used to obtain the ideal mapping of an input subband 
S^i^TT, ^n) on the output band 5*^(0, tt) or Sy{-T:,0) H. Let 

x„(z) ^ [xiz),x{zWQ),--- ,x{zW§-^)f 

ym{z) = [y{z),yizWp),--- ,y{zWj^-')f (3) 



then it can be shown that 



yp(^) = pr(zi/^)Wpy„(zi/^) 

where [r(z)]ij = z'^Sij and [WatJ^j = W^. Substituting from (HI into Q, we get 



(4) 



y„(zi/^) = H„(zi/^)x™(zi/«) (5) 
where K = lcm{P, Q), and Hm(z^/-^) will be given by 



(6) 



H„j(z) is called alias component matrix, or modulation matrix. It can be seen that Q gives a relationship between 
input subbands and output subbands. By proper choice of the dimensions (P, Q) of the modulation matrix H,„(z) 
and by proper selection of the elements of Hm(z) from {0, 1} any required mapping can be obtained. It is shown 
in in that mapping <S'a:(^, ^) ^ Sy{0, tt) can be obtained by using a {P, Q) shift-invariant system, or a P x Q 
modulation matrix, where Q = 2lcm{qi, and P ~ Q — A modulation matrix Hm(z) consisting of 
elements from {0, 1} will be referred as Ideal Modulation Matrix (IMM). 
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In this paper we consider the problem of obtaining the paraunitary polyphase matrices of the filter banks that 
can be realized by {P, Q) shift-invariant systems. In the next subsection we define the problem more precisely and 
highlight the contributions of the paper. 




}. It is required to split the input frequency 
'f-TT, Y^^=^ — tt]. Our objective is to design a 



paraunitary polyphase matrix of minimum dimensions that achieves this splitting. In order to achieve this objective, 
we address the following issues: 

1) If an input subband Sx{^^^^ is to be mapped to output band Sy{0, tt), or Sy{—TT, 0), what would be the 
smallest dimensions of the ideal modulation matrix ? We show in section that though the above mapping 
can be obtained by a P x Q ideal modulation matrix, where Q ~ 2lcm{qi, ^2) and P = Q — as 
given in 13], these are not the smallest dimensions in all the cases. In the same section we obtain smallest 
dimension ideal modulation modulation matrices for a given mapping. 

2) We prove that a (P, Q) shift-invariant system can be characterized by gcd{P, Q) filters. This result as such is 
not new. Many structures that realize {P,Q) shift-invariant systems using gcd{P,Q) filters are given in ||2l- 
However we could not find a method to use these structures to design an RFB based on these structures. In 
section Hn] we obtain a set of gcd{P,Q) filters that is different from the filters given in 121. We show that 
using the ideal modulation matrices, we can obtain ideal frequency spectrums of these filters. These ideal 
frequency spectrums are then used to form an objective function. By minimizing this objective function we 
obtain the required paraunitary polyphase matrix. 



The notations that we use in this paper are as follows. The frequency spectrum of a signal x{n) between oji 
and UJ2 will be denoted by Sx{cdi,oj2)- All vectors will be represented by boldface letters and matrices will be 
represented by upper boldface letters. For two integers k and m, k modulo m is denoted by (fc)„,. For a matrix 
H, the /m-th element is denoted by [H];™. A set of integers from Ni to N2 will be denoted by {Ni..N2}. For a 
polynomial matrix D(z), 'D{z) = D^(z^^), where the * subscript indicates taking complex conjugate of polynomial 
coefficients. x('^i,'^2) denote a frequency spectrum that is 1 for [a;i,W2] and zero elsewhere. 

II. Structure of the Ideal Modulation Matrices 
We recall that the modulation domain equation is given by (|5]l. Putting z — exp{juj) in this equation, we get 



It can be noted that for uj £ [— 7r,7r], the elements of ym(e-''^^™^) cover disjoint segments of tt, tt) and the 
union of these segments is equal to S'j,(— tt, tt). Consider a subset {uji,uj2) of [— tt, tt]. For this subset we define 
an ideal modulation matrix H/(wi,a;2) to be a matrix that maps some segments of Sx on the respective segments 
of Sy according to given desired mapping when oj G {oji,lu2)- We say that a mapping is possible if we can obtain 
disjoint subsets of [— 7r,7r] such that their union is [— 7r,7r] and for each subset, we can obtain an ideal modulation 
matrix. 

Example : Consider the following mapping 



C. Notations 



ym(e- 



(7) 
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Sy{0,n) 



^ Sy{-7r,0) (8) 
For this mapping we have H/(— tt, 0) and H7(0, tt) are given by 

" 1 
H/(-7r,0) = 10 
1 

" 1 

H/(0,7r) = 1 (9) 
1 

It can be seen that y™(eJ"/3) = H/(-7r, 0)x™(eJ"/5) ^^lem w G (-7r,0) and ym(e^"/3) ^ H7(0, 7r)x™(e^"/5) 
where uj S (0, tt) together give the required mapping. 

We now consider the following problem - given two rational numbers ^ and where gcd{pi,qi) = 1 and 
gcd{p2,q2) = 1, find out the the smallest modulation matrices that achieve one of the following mappings: 



'-'x( , ) — ? 

91 92 



Sx{- 



92 



91 



or 



PlTT P2TT 

91 92 



Sxi- 



P2Tr PlTT 



92 



91 



5,(0, tt) 
Sy{-7r,0) 

Sy {-7r,0) 
5,(0, tt) 



(10) 



(11) 



Let Q = lcm{qi,q2). Pi ^ ^ and P2 = Let P = P2 - Pi, m = gcd{P, Q), p = P/m, q = Q/m. The 
following lemma gives the smallest ideal modulation matrices. 

Lemma 1: The ideal modulation matrices will have dimensions P x Q if and only if either of the Pi or P2 is 
even. If both Pi and P2 are odd, then the ideal modulation matrices will be of dimensions 2P x 2Q. Moreover the 
ideal modulation matrices are given as follows- 

1) Pi is even: In this case we will have two modulation matrices H/(— tt, 0) and H/(0, tt). These matrices are 
given as follows - for I ^ [0, P — 1], [H/(— tt, 0)];fe = 1 if fc is given by 

' ^ ^ ' 21 <P 



k = 



(Q-i^ + P-l)) 21 >P 



(12) 



else [H/(— TT, 0)]/A; = 0. H/(0,7r) is obtained from [H/(— tt, 0)]//c as follows 



[H,(0,^)]z, = [HK-7r,0)](p_i),(Q_,)Q (13) 

2) P2 is even: In this case also we will have two modulation matrices H/(— tt, 0) and H/(0, tt). These matrices 
are given as follows - 

for / e [0,P- 1], [H/(-7r,0)];fc = 1 if A; is given by 



(Q-(f-O) 2/<P 



Pi-P 
2 



21 > P 



(14) 
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else [H/(— TT, 0)];fc = 0. H7(— 7r,0) is obtained by using ( fTST l. 
3) Both Pi and P2 are odd: Since the dimension of the ideal modulation matrices becomes 2P x 2Q, this case 
can be considered similar to the above two cases, with Pi, P2, P,Q replaced by 2Pi, 2P2, 2P, 2Q in (fTZt and 
(IHl. 

Proof: We first prove that the number of columns in the modulation matrix must be a multiple of Q. Let 
the dimension of the modulation matrix be i? x 5* and L = lcm{R,S). The modulation equation will then be 
ym{z^^^) = ilm{z'^^^)^m{z^^^)- We first prove that S needs to be a multiple of lcm{qi, q2), i.e. Q. The s- 
th term in Xm(z^/^) and r-th term in ymiz^^^) will be x{z^^^W§) and j/(z^/^W]^) respectively. The term 
x{z^^^W^) covers the input spectrum S^{^^^, e-27^£-i)) f^j. ^ g ^^je ^ ^j{e+27T)y jf [H{z^/^)]rs = 1, then this 
input spectrum maps to the output spectrum Sy{ ^~^^ , y Consider now the mapping given in ( fTOb . In 

order to map the input frequency to the output frequency 0, we require 9,r and s such that 

9_-2n^ ^ Pin ^^^^ 
S qi 

and 

(.« 

Substituting 6 from ( fTSI l into (fT6] l, we obtain + 2(s — r) = 0. Since gcd{pi,qi) = 1, this condition will be 
fulfilled only if 5* is a multiple of qi. By a similar argument it can be shown that S needs to be a multiple of q2 
in order to map the input frequency to the output frequency tt. Thus S needs to be a multiple of both qi and 
q2, or in other words, S needs to be a multiple of lan{qi,q2), i.e. Q. 



If S" = Q, then the condition ^ + 2(s - r) = reduces to Pi + 2(s - r) = 0. Thus if 5 = Q, and the mapping 



Sp 

is given by ( fTOb then we require Pi to be even. By a similar argument, it can be shown that if 5* = Q, and the 
mapping is given by ( fTTT i then we require P2 to be even. If both Pi and P2 are odd, then it is clear from this 
condition that S is required to be equal to 2Q. 

We now obtain H/(— tt, 0) when Pi is even. Consider the frequency in Sy. For 21 < P, this frequency 
will be in 5'j,(— 7r,0). The frequency in Sx that should map to this frequency is 

(17) 



Q Q Q ' 2 
Let k be defined as 

fc = Y + / (18) 

Then if [H,(-^, 0)],^ = 1, it would map ^,(- 12^, _2^) to Syi-^21+P^, -^). 

For 21 > P, the frequency in Sy will be equal to 27r — ^ = 7r(2P-2i) ^ r^j^^ corresponding 5^; frequency 
is given by 

which gives fc as 

k=[Q-{^+P-l)] (20) 
It can be seen that H/(— tt, 0) and H/(0, tt) together achieve the mapping given in ( fTOb . 

The structure of the modulation matrices for rest of the two cases can be obtained in the similar manner ■ 
Remark: Two results close to this lemma appear in H], proposition 3.2, and in [S), theorem 1. However the result 
of this lemma is more general than the above two results and also unifies these results. Proposition 3.2 of H] 
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considers only the case when Pi is even and P > 1. Theorem 1 of [31 always obtains ideal modulation matrices 
of dimensions 2P x 2Q, that are not the smallest ideal modulation matrices if Pi or P2 is even. 

Example 1: Consider again the mapping given in (O. It can be seen easily that since in this case Pi is even, the 
ideal modulation matrices given in (|9]l can be obtained by ( fT2] i. 

Example 2: Consider the following mapping 



3' 



(21) 



Since both Pi and P2 are odd, the ideal modulation matrix dimensions will be 4 x 6, and the ideal modulation 
matrices can be given either from ( fT2] i or from ( fT4] i. Considering (fl4] i for ideal modulation matrices, we obtain 



H/(-7r,0) = 



1 









1 



1 



H/(0,^) 



1 

1 

1 

1 



(22) 



If P is even (i.e. both Pi and P2 are even, or both Pi and P2 are odd) then both the kinds of mappings given in 
([Tol l and ( fTTl i are possible. We now obtain the relationship of ideal modulation matrices for one mapping with those 
of the other mapping. Let the dimensions of ideal modulation matrices he. R x S. Then the modulation equation 
will be 



where K — lcm{R,S). Let P(fc) be the permutation matrix that rotates Xj„(e' 
x^(z) = P(fc)x„(z) then 



(23) 

upward by k. Thus if 



The following lemma gives the relationship of the ideal modulation matrices that give mapping (fTol i with those 
that give mapping (fTTT i. 

Lemma 2: If P is even then an ideal modulation matrix Hj(wi,W2) corresponding to mapping (fTTl i can be 
obtained by the ideal modulation matrix H}(wi,W2) corresponding to mapping (fTOl i as follows, 

where R x S is the dimension of ideal modulation matrices. 
Proof: Consider the modulation domain equation 

y,„(e^"/^)=H}(c.i,c.2)x,„(e^"/^) 

changing oj/R to tt + lj /R, we obtain 

y,„(e^(-+"/«)) = H}(c.i,c.2)P(P/2)x„(e^"/^) (26) 

Let Hj(a;i,cj2) = H}(cji, aj2)P(i?/2). It can be seen that if H}(cji,cj2) gives mapping corresponding to dTol l. 
then H|(cji,cj2) gives mapping corresponding to (fTTT i. ■ 



(24) 



(25) 
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In the following section, we assume without loss of generality that either Pi or P2 is even, since if both of them 
are odd then we can double the dimensions of the shift-invariant system as shown in the above lemma. 



III. Filters Characterizing the Modulation Matrix 
We recall from (|6]) that the modulation matrix is given by 

H™(z^/^) - ^Wt,rp(z-i/^)Hp(z)r„,(zi/Q)WQ (27) 

where K — lcm{P, Q). 

Let m = gcd{P,Q), p = P/m and q ~ Q/m. Then K = mpq. Let (a, 5) be a pair of integers such that 
ap + bq — 1. An element of limiz^^^) is given by 

mp~l 

[Uraiz'/'^nh = — E z'^P'-'^''^/^PmXT"'"Hrs{z) (28) 

"^'^ r=0 s=0 

We now have the following lemma. 

Lemma 3: The modulation matrix can be characterized by m filters. In perticular we assert the following 

1) For a given integer d G {0..m — 1} there are mpq pairs {I, k) such that (/ — k)m ~ d, where I G {O.jjip— 1} 
and k G {0..mq — 1}. 

2) If {I, k) and k') are two pairs of indices such that {I — fc),„ = {V — k'),n then 

[H„(;^i/™P«)]pfc, = [H„(zi/'"P«iy„%)],fe (29) 

for some integer g such that < g < mpq — 1. 

3) g = iff {l\k') = {Ik). 

Proof: In order to prove the first part, for a given I we count all the k that give (/ — k)m = d. Such k are 
given by A: = {{I — c?) + rm)mq, Q < r < q — 1. Thus for a given ^ we have q values of k that give {I — k)m = d. 
Since / takes values from to mp — 1, total number of pairs {I, k) satisfying {I — k)m will be mpq. 
In order to prove the second part, we note that 

mp—lmq — l 

[H™(zi/"f«)],,,, = — E E z'-^W^^l^''^P'-^''-'^'i-WX'-''i-Hrs{z) (30) 

"^'^ r=0 .s=0 

Let [V — l)mp = t and (fc' — k)„iq = t + hm, for some integers t and h, then it can be shown that 

Txr{k'-k)ps-{l'-l)qr „ Wff(ps-'3'-) /tin 
mpq **7npq ^ J 

where <? = (i + hamp). Using this in ( l30b we get 



This proves the second part. 

To prove the third part, we note that if (I' , k') = (I, k) then t = 0, /i = and we have g = 0. Now if g ~ 0, 
then t + ahmp ~ fmpq for some integer /. This gives t = { fq — ah)mp. Thus i is a multiple of mp. But since 
< t < mp — 1, the only possibility is t = 0, and this gives I' = I. In order to prove fc' = fc, we note that with 
t = we have ahmp = fmpq. Now since a is not a multiple of q, h should be a multiple of q. Let h = vq, then 
(fc' — k)mq = vmq, giving v = 0, and fc' = fc. 
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To show that the modulation matrix can be characterized by m filters, for a d G {0..7ti — 1} choose any index 
{l,k) such that {I - k),n = d. Define = [H„(zi/"P9)]/fc. Then from the above results, mpq elements 

of the modulation matrix, namely [H„i(z)]/'fc/ such that — k')m = d can be written as 

[H^(zi/™P«)]z-fc' = i7d(2'/"^W„%) (33) 

In this manner all the elements of H„i(z^/™p'') can be written in the form of i7rf(z^/™P''W^jpq) for some d G 
{0..m — 1} and some g G {0..mpq — 1}. ■ 

By using the lemma([T]i and the lemma© we can obtain ideal frequency spectrum of the m filters that characterize 
the modulation matrix. The following lemma gives ideal frequency spectrum of these m filters: 

Lemma 4: Only maximum two filters out of m filters will have non-zero ideal frequency spectrum. Moreover 
by appropriate choice of the filters, it can be shown that the filters having non-zero ideal frequency spectrum will 
have frequency spectrum given by x(;^, ■^^^^^) and x(^^^i^p^, ~;^) for ^0™^ integer v. 

Proof: In order to prove that only two filters out of m filters have non-zero frequency spectrum, we note that 
when Pi is even, then it can be seen that if H/(— 7r,0)];fe = 1, then either [l — fc)„i = (Pi/2),„ or [l — fc),„ = 
(-Pi/2),„. Let Hdi{z) and Hd^iz) be two filters out of m filters such that [H„i(zi/™P«)],,fc ^ iJdi(zi/"P«W,9p^) 
when (/ - fc)™ = (-Pi/2),„ and [H„,(z1/™p<z)],_,. ^ i/d2(^^/™P'M/,?,p,) when {I - fc)„ = (Pi/2)„„ then it can 
be seen that the ideal frequency spectrum of only Hdi{z) and Hd2{z) will be non-zero. 

Similarly if P2 is even, then Hdi{z) and Hd2{z) will have nonzero ideal frequency spectrum if [Hm(z^/™''*)];,fe = 
Hdi{z^'"'P''Wlp^) when (/ - fc)„ = (-P2/2)™ and [H™(z1/™p9)];^^ = iJ^aC^^/^PW^p,) when {I - fc)™ = 
(P2/2),„. Rest of the filters will have zero ideal frequency spectrum. 

In order to prove the second part, we consider two cases 

1) Pi is even; For this case if we choose Hdi{z) and Hd2{z) such that 

Hdi{z^/^^'>Wl^^) = [H™(^'/'"'"^)]o,p,/2 (34) 
i/.2(z^/"^^W^-;,) = [H,„(^^/™^')]o,(™,_i^) (35) 

then it can be shown that 

Q I \ -{2v + P)tt ~2im 

<SMi(-7r,7r) = x( , ) 

mpq mpq 

S^hd2(-7r,7r) = x( , ) (36) 

mpq mpq 

If {Pi)m = then we have dl = d2 and for such a case 

mpq mpq mpq mpq 

with V = 

2) P2 is even: For this case if we choose Hdi{z) and Hd2{z) such that 

iJ,2(^^/"^«W^-;,) = [H,„(zi/™^^)]o_(„,_i^) (39) 

then it can be shown that 



(37) 



(38) 
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~2vn -{2v_P)7r 



mpq 



mpq 



A2v-P)tt 2vk . 



mpq mpq 

Again if {P2)m = then we have dl = d2 and for such a case 



<5'/idl( — TT, Tt) 



, -2vt: -(2u-P)7r, ,(2v~P)7r 2v'k , 



mpq 



mpq 



mpq 



mpq 



(40) 



(41) 



with = 



The ideal frequency spectra of these filters are used to form an objective function, minimizing which would give 
us desired paraunitary polyphase matrix. 



IV. Design of the Paraunitary Polyphase Matrix 

A. Formation of an objective function 

Given rational numbers ■ • ■ > }, our objective is to design a paraunitary polyphase matrix Hp(z) that 
gives the required splitting of the input frequency spectrum S'^(0, 7r)into bands. 

Consider n-th channel of such a filter bank. The part of the input spectrum corresponding to this band is 

It is clear from lemma[T]that the ideal modulation matrices for this channel will be of the size P„ x Qn or 2P„ x 2Qn- 
Let the dimensions of ideal modulation matrix be denoted by x Qn- 

Since we have to realize a paraunitary polyphase matrix for the filter bank, it is required that the polyphase 
matrices for all the channels should have same number of columns. Let S = lcm{QQ,Qi, ■ ■ ■ ,(5Ar-i),i?n = 
PnS/qn,Kn = J^jt^o ~ '^]=o ~ 9cd{Kn, S). Then the polyphase matrix for n-th channel, denoted 

by Hp„(z), will be of the dimensions Kn x S and the polyphase matrix for the filter bank will be given by 



Hp(z) 



Hpo(z) 
Hpi(z) 



H 



.p(JV-l) 



(^) 



(42) 



Let the modulation matrix corresponding to Hp„(z) be denoted by Hto„(z). The modulation matrix H„i„(z) 
will be characterized by M„ filters. Let these filters be denoted by Hin{z),0 < i < Mn ~ 1- Since Hp„(z) is 
paraunitary, from ^ it can be seen that Hm„(z) is also paraunitary. By using the paraunitary condition on the 
modulation matrix Hm(z), i.e. Hm(z)Hm(z) ~ I, it can be shown that the filters characterizing the modulation 
matrix satisfy 

M„ - 1 

|2 



E 

1=0 



KUen\ =Kn 



(43) 



where Kn is some constant. From this, it can be seen that if the stopband energy of all the filters is minimized, then 
the passband energy would get maximized. Thus we can define an objective function, based on stopband energies 
of filters in all the channels. This is done as follows: 

1) Find out the dimensions of the polyphase matrix Hp(z) as explained above. 
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2) Define a paraunitary polyphase matrix Hp(z, Q) that depends on the parameters Q. 

3) Find out the ideal frequency response for all the filters Hin{z), < i < A/„ — l,0<?7<iV— 1. 

4) Choose a transition bandwidth e. Using this transition bandwidth and ideal frequency response of iJm(z), 
find out its stopband 7™ 

5) Define the objective function 

\H^n{en\^ dio (44) 

,1=0 1=0 -^T*" 

By minimizing D{Q) with respect to 8 we can obtain the required polyphase matrix tip{z, Q). 
B. Design Examples 

For the design examples given below, we will use Given's factorization for the paraunitary matrices (|5|. Thus 
if a paraunitary matrix has size N x N and order K, then the parameter vector 6 will have K{N — 1) + NC2 
parameters. 

Example 1: In this example we consider the rational splitting {|, i, |}. 

Since for the third channel we have Pi = 3 and P2 = 5, accordng to lemma(|T]l the size of polyphase matrix for 
this channel will be 4 x 10. Thus the complete filter bank can be realized by a 10 x 10 polyphase matrix, with first 
4 rows for first channel, next 2 rows for second channel and last 4 rows for the last channel. 

By applying lemmaQ and lemma© it can be seen that each channel can be characterized by two filters and 
only one filter out of these two filter has non-zero ideal frequency spectrum. 

A polyphase matrix is realized using the procedure outlined above with transition bandwidth equal to pi/ 20 
and the number of stages in the polyphase matrix equal to 7. The ideal frequency spectrum and actual frequency 
spectrum of the filters are given in Fig.©. In these figures stopband of ideal filters is shown at —60dB line. 

Example 2: In this example we consider the rational spUtting {|i 5, 

The size of the polyphase matrix in this case would be 9. First channel will have only one filter, next two channels 
will have three filters and last channel will have only one filter For second and third channel, two filters will have 
nonzero frequency spectrum. The ideal and realized filters are shown in Fig.©. For this realization the transition 
bandwidth is taken as 7r/20 and the polyphase matrix has 7 stages. 

V. Conclusion 

In this paper we have presented a design method to design rational filter banks using {P, Q) shift-invariant 
systems. We have obtained minimum dimension of a (P, Q)-shift invariant system for ideal mapping. Next we have 
obtained a set of gcd{P, Q) filters characterizing (P, Q) shift-invariant systems. Ideal frequency spectrum of these 
filters are obtianed obtaining and then paraunitary polyphase matrix of the filter bank is obtained by minimizing 
the stopband energies of these filters. 

Since the design method is based on (P, (5)-shift-invariant systems, filter banks with arbitrary rational splitting 
can be designed by this method. Further improvements in the method are possible by finding out a good "starting 
point" for objective function minimization. 
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Fig. 3. Filters for three channels of the filterbank in Exanple 1 - Ideal frequency spectra are also shown with actual frequency spactra. The 
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[3] R. G. Shenoy. Multirate Specifications via Ahas-Component Matrices. IEEE Trans on Ckt. and Sys.-II, 45(3):3 14-320, 1998. 
[4] E. W. Bai T. Chan, L. Qiu. General Multirate Building Structures with AppHcation to Nonuniform Filter Banks. IEEE Trans on Ckt. and 
Sys.-II, 45(8):948-958, 1998. 

[5] M. Vetterli and J. Kovacevic. Wavelets and Siibband Coding. Prentice Hall PTR, Englewood Cliffs, New Jersey, 1995. 



11 
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Fig. 4. Filters for Example Two - Ideal frequency spectra are also shown with actual frequency spectra. The stopband of ideal frequency 
spectra is shown at -60dB 
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